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Abstract. We present a geometric approach to the three-body problem in the 
non-relativistic context of the Barbour-Bertotti theories. The Riemannian metric 
characterizing the dynamics is analyzed in detail in terms of the relative separations. 
Consequences of a conformal symmetry are exploited and the sectional curvatures of 
geometrically preferred surfaces are computed. The geodesic motions are integrated. 
Line configurations, which lead to curvature singularities for N 3, are investigated. 
None of the independent scalars formed from the metric and curvature tensor diverges 
there. 



1. Introduction 

The nonrelativistic dynamical models of Barbour and Bertotti || arose from the 
criticism of the concepts of absolute space and time. They describe a classical 
interacting TV— particle system subjected to the Hamiltonian, momenta and angular 
momenta constraints. The invariance group of the theory is the Leibniz group 
[01, which includes time-dependent translations and rotations together with the 
monotonous but otherwise arbitrary redefinition of time. 

In a previous paper 1^ being referred hereafter as paper I, one of the present 
authors has analyzed the underlying geometry of the Barbour-Bertotti theories. The 
reduction process on the Lagrangian was carried out by solving the constraints, 
arriving to a Riemannian line element. This reduction was possible for all 
configurations but the line ones. The emerging Riemannian metric G was shown 
to represent the first fundamental form of the orbit space of the Leibniz group. The 
geodesies in this metric characterize the free motions, pertinent to constant potential 
V. For a generic potential V ^ const the motions are geodesies of the conformally 
scaled Jacobi metric —2V{x)G. 

In I the Riemann tensor and curvature scalar were computed in terms of the 
vorticity tensors of the generators of rotations. Then the curvature scalar was 
expressed in terms of the principal moments of inertia and the number of particles. An 
analysis based on this expression allowed us to conclude that the line configurations 
represent curvature singularities for iV 7^ 3. One would like to say more about these 
configurations for the exceptional case N = 3. This is one of the motivations of the 
present work. 
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The second motivation for specializing to three particles is the remarkable 
coincidence between the number of relative separations N(N—l)/2 among the particles 
and the dimension 3N — 6 of the reduced space | . This feature enables one to employ 
the distances as a symmetric set of variables of the space of orbits and to analyze in 
detail the underlying Rieniannian geometry. 

We discuss and picture the space of orbits in Sec. 2. Similar discussions can be 
found in g and §]. The space of orbits being a manifold with boundary, we announce 
and prove the conditions a geodesic reaching this boundary has to obey. 

The rest of the paper is organized as follows. First we particularize the generic 
expression of the curvature scalar obtained in I to the case of three particles. For this 
purpose we compute in Sec. 3 the explicit expressions of the principal moments of 
inertia in terms of the relative separations. We find that the curvature scalar reduces 
to a particularly simple form in the case of three particles. 

We reveal more details of the geometry in Sec. 4. There we give the reduced metric 
on the space of relative separations first in terms of the distances and second in terms 
of a radial coordinate and suitably chosen angular variables. Using the second set of 
variables the Riemann and Ricci tensors are computed. None of the scalars formed 
from the metric and the Riemann tensor diverges in the line configurations. 

We show that the metric is conformally flat. Furthermore, the metric has a 
conformal symmetry. This enables us to compute the extrinsic curvature of the 
ellipsoid surfaces orthogonal to the conformal Killing vector in Sec. 5 and to 
demonstrate in Sec. 6 that the radial lines are geodesies. 

The purpose of Sec. 5 is to analyze the eigenvalue problem for the Ricci tensor. 
Remarkably, one of the Ricci principal directions is the conformal Killing vector. 
The sectional curvature of the ellipsoid surfaces orthogonal to the conformal Killing 
vector is found by the Gauss relation. The sectional curvature of the conical surfaces 
orthogonal to the other eigenvectors are also computed. All eigenvalues of the Ricci 
tensor are finite for the line configurations. 

Free motions are shown in Sec. 6. to be geodesies of the reduced space. The 
geodesic equation is integrated in the space of distances. 

Throughout the paper we use the following notations. Latin indices denote 
components in the space of distances. They run from 1 to 3 and are raised or lowered 
with a non-flat metric. Greek indices label quantities pertinent to different eigenvectors 
of the Ricci tensor. There are few exceptions under these rules, indicated where 
necessary in the text. Summations are explicitly written whenever the summation 
convention can not be applied. Partial derivatives are denoted by comma. 

2. The space of orbits 

The nine coordinates characterizing the positions of the particles can be chosen as the 
relative separations Ui = (a, b, c), the Euler angles = {a, f3, 7) of the normal to the 
plane of the three particles and the coordinates of the center of mass x\ By freezing 
the translational and the rotational degrees of freedom, we fix the coordinates of the 
center of mass and the Euler angles respectively, such that the distances (a, b, c) will 
coordinatize the reduced space. Before proceeding with the analysis of the metric 
and curvature properties of the reduced configuration space, we briefly describe and 



f 3A'' — 6 = N(N — l)/2 holds also for A'' = 4. A discussion of the adequate coordinates for this case 
can be found in M. 
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Figure 1. The intersection of an r = constant surface with the space of orbits for 
three particles in the Cartesian frame of the distances is a spherical triangle, the 
boundaries of which are arcs of great circles. 

The domain of the polar and azimuthal angles 9 and (p defined in the Cartesian 
system spanned by the coordinate lines is restricted first by the positivity of af. 

< < 1 , < ^ < I , (1) 

second by the triangle inequality: 

1 sin (f — cos (p I < cot 6 < sin cp + cos ip . (2) 

Therefore the reduced space can be imagined like a diverging beam of rays starting 
from the origin, with a triangular section, as can be seen on Fig 1. The boundaries 
of the admissible domain of the reduced space are defined by the equalities in Eq. (||) 
as planes passing through the origin. In consequence each constant radius section of 
the reduced space is bounded by arcs of great circles. The origin O corresponds to 
the unphysical situation of all three particles in the same location while the planar 
boundaries to the situation of the three particles on a line. The edges of the beam 
(the intersection of two plane boundaries) are again unphysical: there the positions of 
two of the three particles coincide. On all boundaries of the space of orbits the inertia 
matrix has vanishing determinant, consequently the metric (32) of I is ill-defined there. 
On the other hand the scalar curvature is well behaved in the line configurations, 
suggesting that the boundary planes (including the edges but not the origin O) are 
merely coordinate singularities. Topologically the region inside the boundaries is the 
quotient space R^/E{3), where E{3) — xi SO{3) and is homeomorphic to half of 
R^. The boundary is homeomorphic to R^ Q. 

The special situation of the restricted three-body problem, with a lighter third 
body orbiting about the other two heavier, whose separation is not changing is 
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confined within intersections of tlie space of orbits (the "beam") with one of the 
planes = const. 

As physically the line configurations are no more special then others, our 
description has to allow the geodesies to reach and to leave the boundaries. We then 
require that a geodesic reaching the boundary at some instant to should not leave the 
admissible domain of the {ai} space, by imposing the following boundary conditions: 

a±(to) = , (3) 
a±(to)<0. (4) 

Here aj^ and aj^ represent the projections perpendicular to the boundary of the first 
and second derivatives of the position vector a — (a, b, c) with respect to the geodesic 
parameter. The perpendicularity is meant in an Euclidean sense, by projecting a and 
a to the outer normal to the boundary surface. 
Proof of the boundary conditions: 

It is well known that the motion of a three body system with vanishing angular 
momentum L — ^j^T-iA = f is confined to a plane By using Cartesian 

coordinates {x, y, z) we choose za = Q = za and then the constraints "YIia ^ Az = 

l^Ax = '^Ay = are trivially satisfied. The distances ai are given as: 

= {x2 - x^f + {y2 - ys)'^ , 
b"" ^ {x3 - xif + {y3 - yif , 

c' = {xi - X2f + {yi - y2f ■ (5) 
As parameter of the geodesic the time can be chosen. The time derivative of (5) gives 

aa = {x2 - X3){x2 - X3) + {y2 - y3){y2 - ys) , 
bi} = (2:3 - xi){x3 - ±1) + {ys - yi){y3 - Vi) , 

cc = (xi - X2){xi - X2) + (yi - y2)(yi - 2/2) ■ (6) 

From among the "velocities" xi, X2, is, yi, 2/2, ys only three are independent due 
to the remaining constraints '^a'Pax = X^A-f'^a ~ '^a^Az = 0. However for the 
purpose of the proof it is not necessary to make this manifest. We just comment 
that in general {ctj} depend on three arbitrary velocities and in the majority of cases 
they are independent. A notable exception is given by the coUinear configurations, as 
will be seen later. In the computation of the second derivatives we employ that the 
motions are free, xa ~ jjA — 0: 

aa + = {±2 ~ ±3)^ + (y2 - ys)^ , 

bb + b^ = {±3 - iif + im - , 

cc + d^ ^ (±1 - ±2^ + (yi - y2 f ■ (7) 

Then we study the coUinear configurations. By a rotation in the xy plane, it can be 
achieved that the coUinearity arises on the a;-axis. We label the particles as (3, 1, 2) in 
order of increasing x-coordinate. In this way X3~ xi = —b, x\ — X2 — ~c, a;2 — X3 — a 

t Here and later in Eqs. ( |ll[ ) and ( p^ ) we use the conventions of I: capital Latin indices count the 
particles and lower case Latin indices their coordinates. Summation is explicitely indicated only over 
particles. 



and the boundary is expressed hy a — b + c, with the outer normal n=(l, —1, —1). 
The system (6) becomes 

d^{x2-X3), b^-{x3-xi), C^~{xi-X2)- (8) 

It is immediate to show aj_ = a • n — 0. Employing Eqs. (|^), the system (7) can be 
written in the form 

(6 + c)d = (2)2 - 2)3)^ , bb = (?/3 - yif , cc = {yi - 2)2)^ ■ (9) 

By simple algebra we find then 

.. " „ Hvi -2)2) + c(j/3 -yi)? ^ ^ 

aj^ = a • n = — < (10) 

bc{b + c) ^ ' 

which completes the proof of the boundary conditions, Eqs. (3)-(4). 
Q.E.D. 

3. The curvature scalar for three particles 

In this section we work out the expression for the curvature scalar (given in I in terms 
of the principal moments of inertia) for the particular case of three particles. For 
this purpose first we write the principal moments of inertia in terms of the relative 
separations. 

By choosing the positions of the three particles in some initial coordinate system 
(Fig. 2) originating in the center of mass as a; = (a;, 0,0), = (^,y,0) and 

o o 

x^'^ = {~X3, —y3,0), the X3 and 2/3 coordinates are determined by the condition 

o 

X3 = 2/3 = — y ■ (11) 

7713 "73 

The remaining three coordinates (x,^, y) are related to the relative distances as: 
TOgfl^ = [mix + (r?72 + ?7i3)^]^ + (r?72 + mTify^ 

77136 = [(rni + 7773)2; + 7772^] + 7772?/ 
2 _ ^2 I „,2 I e2 



C = X 



y^+e^'^xi. (12) 



We introduce the positive quantities A, W and Z , depending only on relative 
distances and with the dimensions of moments of inertia, as follows: 

777i?772C^ + •m2m^a^ + m^mib^ 7712W3 ? 

W — — — , A = a , 

2M 7772 + 7773 

^2 ^ (a + b + c)(a + b - c)(a - b + c){~a + b + c) . (13) 

4iV7 

Heron's formula relates the quantity Z to the area A of the triangle: Z = 
2y^^77ii7ri27n3/Af . From Eqs. (12) and Eqs. (13) a useful relation is found: 

ymAX^'x^''5cd = '2W (14) 

^ ^ o o 



^2(^ >y>o) 




(x, , 0) 



mj(-Xg, -y^, 0) 



Figure 2. The three particles in the center of mass frame. 



A computation employing the definition of the inertia tensor in the center of mass 
frame (Eq. (21) of I) and Eqs. (12), (13) and ( p^ ) yields the expression of the tensor 
of inertia in terms of distances: 



hi 



= 1 -^hi{2W 
.0 



hi) 



-^hi{2W 

Z2 2W - hi 




hi) - 



where 



hi 



{2W - A) 

From here we find the principal moments of inertia 

/3 



h.2 = W± \/W^ - Z2 , 



2W . 



(15) 



(16) 



(17) 



They contain the relative distances and the masses in a symmetric fashion {hi, which 
is the only term containing the non-symmetric expression A, dropped out from their 
expression). By insertion of the principal moments of inertia in Eq. (43) of I we have 
an independent check of Eq. (^4|) . 

The curvature scalar of the reduced space given by Eq. (51) of I, written in terms 
of the distances takes the remarkably simple expression: 

As expected from the argumentation in I, the curvature scalar is well behaved even in 
the case of all three particles in a line. This suggests that in the line configurations a 
coordinate singularity occurs. By contrast, the unphysical situation of all particles in 
a point represents a true curvature singularity, as W vanishes there. 
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4. Metric properties 

In order to find the metric in the space of orbits in terms of the distances we start 
from the degenerated metric (32) of I. Then we pass to the basis associated with the 
coordinates (?/*, aj, x*^), where ?/' — {x, ^, y). The relevant part of the transformation 
matrix is given by 

The 3x3 block with the basis in the 1-forms dy^ then is transformed to the basis dui 
through the relations (12). Finally we get: 

/ m2 + 7713 + 6^ — — 6^ + \ 



m2m^ 2m^ab 2m2ac 

a? ^b"^ — (? mi + nis —a? + 5^ + 



(20) 



2ra^ab mim^ 2mibc 

a? — b'^ -\- (? —c? + 6^ + TOi + TO2 

\ 2m2ac 2mibc mim2 J 

The covariant form of the metric has a more complicated expression. Eq. ( pO| ) 
is the induced metric in the space of orbits. It is not surprising, that G*-' is 
manifestly symmetric in both masses and distances, because the overall translational 
and rotational degrees of freedom have been suppressed. 

Remarkably, the metric G^^ has a conformal symmetry. Indeed, ( pO| ) is unchanged 
under multiplying all distances with the same factor, thus G"^ is homogeneous of degree 
zero in a^. The Euler theorem for homogeneous functions gives then G'^'^.a^ = 
for any point of the reduced space with coordinates a^. Then the vector 

K' = a, = (a, b, c) (21) 

is a conformal Killing vector, generating a homothctic motion: 

Ik G'^ = -2G*^ (22) 

The norm and covariant components of this conformal Killing vector can be expressed 
in terms of the previously introduced quantity W: 

K = V2W , Ki = -^(777277730, mimsb, mim2c) = V^W . (23) 

Here V denotes the connection compatible with the metric G. 

The distance in the metric Gij of any point from the origin is the length of the 
corresponding conformal Killing vector (|2^). Thus the curvature scalar (|l^) is just 6 
divided by the square of this distance. 

Whenever one has a metric G'^ with a conformal symmetry, by an appropriate 
rescaling a new metric G*-' = Vl^G^^ can be found which has the property, that the 
conformal Killing vector of G'-' is a true Killing vector of the metric G*-' . In our case 
the scaling function was found to be fi^ = 0? + b"^ + (? — . Choosing an adapted 
coordinate system the metric G'^ does not depend on r. Such a coordinate system is 
given by {In r, 9, ip) in which the conformal Killing vector of G'-' takes the simple 
form if* = (1,0,0). Thus the metric G'^ is expected to have the form: 

= L . G'^angles) . (24) 
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Unfortunately the attempts to write the metric in terms of the radial coordinate r, 
supplemented by any two angular coordinates, have resulted in less symmetric forms. 

In the search for angular variables which give a reasonably simple form of the 
metric, we have found a convenient (but redundant) set given by the sines 

Si — sin 9i (25) 

of the three angles 6i (Fig. 2.). The condition di + 02 + 63 — tt constrains the variables 
Sj as follows: 



sl + sl-sl^2siS2^^1-4 (26) 

and all of its cyclic permutations. The square of Eq. (p6| ) yields a symmetric form of 
this constraint: 

(Sl + S2 + S3)(S1 + S2 - S3)(S1 - S2 + S3)(s2 + S3 - Si) = 4:slslsl . (27) 

From the sine theorem we find that the relative separations are related to the 
redundant variables Si as: 

T 

ai = —Si , (28) 

where S'^ = sl + S2 + S3. By multiplying Eq. ( p7| ) with (r/S)^ we find S'^ in terms of 
the distances: 

^ ^ 91,2 2 ■ (^^) 

We also introduce the quantity 

2 2MWS^ 2 2 2 2M^WZ^ 

a = 5 = mim2S3 + m2m3Si+m3miS2 = „,„ „ , (30) 

r7iiTO2r7i3a^o^c"^ 

which depends only on the angular variables Si as can be seen from the second form 
given in Eq. (pO|). The third form arises by inserting Eq. ( p9| ) in the first form and 
similarly to Eq. ( p9[ ) is useful in rewriting the angular expressions in terms of the 
distances. 

By multiplying Eq. (Eq) with {r/S)"^ the following relation emerges: 



Similar equations follow for the cyclic permutations. Therefore the contravariant 
metric expressed in these new angle variables, while keeping the basis of the metric in 
the 1-forms c?a, , takes the concise form: 



k 

Here jjL^^^ is the "square" of the symbol e*-'*^, zero for any pair of coinciding indices 
and taking the value 1 otherwise and no summation over j applies. 

The concise form ( p^ of the metric allows us to study its signature. First we 
remark that all three diagonal elements of G^^ (the first term in (p2)) are positive. Also 
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all three sub-determinants of the diagonal elements have similar forms, for example 
the sub-determinant of G^^ is 

+ > 33 

and are also positive. Finally the determinant 1/G of G"^ , 

= (34) 

is positive except in the singular configuration of particles in a line, where is zero. 
Therefore we are in the position to conclude that the metric is positive definite, modulo 
the boundaries where it diverges. 

Direct computation of the curvature scalar from the metric yields (^8|) once again, 
which in terms of the angular variables and r takes the form: 

6 A/52 

R = . 35 

The Ricci tensor can be expressed in terms of the curvature scalar: 
TO17712TO3 



(36) 



As in three dimensions the Weyl tensor vanishes, the Riemann tensor is determined 
completely by the Ricci tensor, the curvature scalar and the metric: 

Rijki = '2Gi[kRi]j - 2Gj[kRi]i - RGi[kGi]j . (37) 
Again we write the result of the computation in terms of the curvature scalar: 

Rijkl = R ^^^^ 2^ ^ijaSa ^klbSb ■ (38) 

a b 

All multiplying factors of R in Rijki and depend only on the angular variables, 
thus these tensors depend on r only through R. 

Next we study the behavior of the curvature in the line configurations. Though 
the curvature scalar is undetermined in the form (|35|), it is clearly nonsingular when 
written in terms of the distances (p^, as already remarked. From Eq. ( ^8| ) we compute 
the Kretschmann scalar and find that it is nonsingular either: 

R^JklR''''' = R^ . (39) 

From Eq. (|3^) an algebraic relation can be deduced between these scalars 

R^,klR''''' - ^R^jR'' + i?' = , (40) 

which gives 

2 

In fact in three dimensions there are only three independent scalars which can be 
formed from the Riemann tensor and the metric. These can be given either as 



R,,R!= = — . (41) 
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or as the roots of the secular equation Det{Rij — XGij) — To complete the first 
, (43) 



set (42) we compute 
Det{R) 



Det{G) 

while the secular equation will be dealt with in the next section. 

Finally we remark that in three dimensions the vanishing of the tensor 

Rabc = ^cRab - ^bRac + ^{Gac^bR - Gab^cR) (44) 

is equivalent with the conformal flatness of the metric The computation gives 
Rabc = 0, therefore the metric is conformally flat: 

Qab ^ ^2^ab ^45) 

Here rjab denotes a flat metric. The conformal factor can be determined by solving the 
differential equation which arises from the comparison of the expression (|l^) with the 
formula relating the curvature scalars of the metrics G""^ and Ty"** (the curvature 
of the latter being zero) : 



R= ir" - :l ' ■ (46) 



It is easy to prove that 

^'^ ^aW , (47) 

a arbitrary constant, solves the differential equation, by employing Eqs. ( p3[) and the 
relation VaK"' = 3, which stems out from the conformal Killing equation (|22|). 

A coordinate system {wi}, related to the Jacobi coordinates, frequently used |^ 
in molecular dynamics, has the advantageous property that 77°^ = diag{l, 1, 1). The 
price one has to pay in using such coordinates is that they depend on the distances in 
a nonsymmetric fashion: 



wi = 2{B -W) , W2 = 2^/2BW - B'^ - , W3 = 2Z , (48) 

where the notation B = mim^b'^ / {mi + 777,3) was introduced. In the coordinates {wi} 
the constant in the conformal factor ( p7| ) is a = 8. 

5. Sectional curvatures 

In this section wc study the eigenvalue problem of the Ricci tensor (|3^): 

R'^n^=pn'. (49) 

a a a 

The eigenvalues P are and a double root equal to R/2. All eigenvalues are well 

behaved for the configuration of particles in a line. The conformal Killing vector 
K"^ ( [2l| ) is the eigenvector corresponding to the zero eigenvalue. A two-dimensional 
subspace spanned by the vectors: 

A:(\-| (0, c/ms, -6/TO2) 

^(2) = {-c/ms, 0, a/777i) 

^(3) = ib/m2, -a/mi,0) (50) 
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corresponds to the degenerate eigenvalue R/2. Any two of these vectors are hnearly 
independent, but not orthogonal to each other. From the general theory of symmetric 
matrices we know that the two-dimensional space spanned by X'l^^-^ is orthogonal to 
. A compact notation for the eigenvectors is: 

K' = m. ^(V)-E^-^-^ (51) 

and for a dual set of the eigenvectors: 

Ki = ^iijkajmjmk 

^i''^ = ^Yl t^^^j^jarttu + E fi^jkirijal . (52) 

Each of the one-forms f/^'' is the dual of some vector Y^^^, orthogonalf to both K'^ 
and Xl^y 

It is immediate to check that the set of vectors F^'^^j and are surface orthogonal, 
as they satisfy the relations: 

K.^jKke'^^ = , Y^Y^'^^e'''' = . (53) 

We proceed to find these privileged surfaces. It follows from (|3|) that the ellipsoids 

F = W- const = , (54) 

are the surfaces orthogonal to K^. Then the surfaces /^''•' = orthogonal to the vectors 
Y^^-^ can be determined by imposing the proportionality of their gradients with 1^^'^' : 

r/'') = A('')v,/('') , (55) 

by some undetermined factor A^"^. The detailed derivation is given below for the 
surface orthogonal to Y^g^ . We look for this surface in the form 

=e(c)-<i>(a,5) . (56) 

Comparing the expressions of f/^'' from (52) and ( p5|) we find three equations to 
determine the unknown functions Q{c), ^{a,b) and X^^\ which are: 

A(3)$^=m2ac (57) 
A(3)$^2 = mibc , 

By inserting A*^^^ expressed from the first equation into the rest of the system (57), 
the remaining equations can be brought into the form 

c dQ / ,2 2\ d,^ 

-— = a = (777l6 +77720 — jT' (58) 

2 dc d{m\V + m2a'^) 

t The vectors Y^'^j themselves, however have cumbersome expressions, therefore we will avoid to 
write them. They can be obtained by the Gram-Schmidt orthogonalization procedure. The vectors 

Y} , are not orthogonal to each other either. 
(") ° 
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which has the immediate sohition Q — a ln(7c^) and $ = a ln[/3(TOi6^ + m2a^)]. Here 
a, P and 7 are constants. Thus the equation of the surface orthogonal to the vector 
y^g-j is given by: 



= a hi 



7c 



, (59) 



or equivalently as 

= m,ib^ + m2a^ - const • = , (60) 

Similar equations with cyclically permuted quantities define the surfaces to which the 
other two vectors F^'^-j are orthogonal. These surfaces are cones with the tips in the 
origin and they have elliptical sections. 

Next we derive the sectional curvatures of the privileged sections F and F''^^ . We 
pick up an independent set of eigenvectors n * by normalizing and any two of the 

three vectors F^g^ . For each of these surfaces we define the Riemann tensors '^'^^R "^^^ 

like in Eq. (46) of I. We also define the projector operators PjJ and the extrinsic 

a 

curvatures K ah for the sections: 

Pt = 6b-n''nt, Kab^PaPb^cnd ■ (61) 

a a a a a 

The Gauss equation connects the curvature tensor and extrinsic curvature of the 
sections to the curvature tensor of the space of distances: 



^'^i? abed ^ Pi pi Pi PdReSah +2Ka[cKd]b- (62) 

a a a a ^ ^ 



A double contraction with the induced metric on the sections G — n°- n" gives the 

a a 

equations for the sectional curvatures: 

^^^R ^ R-2P +{K\f ~ KabK'''' . (63) 

a a a a a 

Here we have employed that rt ° is eigenvector of the Ricci tensor, Eq. ([49[). 

a — 

Straightforward computation based on the above prescription has shown that the 
sectional curvature of the conical surfaces F^'^^ — vanishes. 

For the ellipsoid surfaces with the normal = \J2W the extrinsic curvature 
is found readily from the equation of the homothetic motion (p^). In terms of the 
covariant form of this equation is: 

Vi^j + Vj-rii = (Gy - n.Uj) . (64) 

Contracting twice with the projectors to the ellipsoid surface, the desired result 
emerges: 

Kii, = -^/= - n^Uj) . (65) 
A simple computation then yields 

{K\f-KabK'''^^^ (66) 
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and from (|63|), the sectional curvature of the eUipsoid surface F = is found: 

(67) 

This is constant on the eUipsoids, therefore in the metric the F = surfaces are 
spheres. 

6. Geodesic motions 

The equation for the homothetic motion ( p^ ) carries even more information. A 
contraction of its covariant form with gives: 

K^VjK, + ]^{K,K')^, = 2K, . (68) 

The second term is simply Ki as can be seen from (p3|). Thus the integral curves of 
the conformal Killing vector field are geodesies: 

K^VjK, = K, . (69) 

We recall that K^d/dx'' = 5/9 In r. Therefore Inr is not an afhne parameter, however 
it is easy to check that pr + q is, if p and q are constants. 

Eq. ( |69| ) says that the motions along radial lines of the reduced space are geodesies. 
Therefore the boundaries of the reduced space (Fig. 1) are geodesic planes in the sense 
that they are spanned by geodesies passing through the origin. The privileged conical 
surfaces F'^"^ = are geodesic surfaces in the same sense. 

The ellipsoid sections F = are not geodesic surfaces. To see this, we compute 
the Gaussian curvature /C of a geodesic surface whose tangent space in each point is 
spanned by any two of the (non-orthogonal) eigenvectors Y^^^ by means of [D 

We have found /C — R/2. Since the Gaussian curvature is half of the curvature 
scalar, the latter would be '^^''Rgeod = 2/C = i? for a geodesic surface with tangents 
Yj-'^j . This is different from the scalar curvature ( |67| ) of the ellipsoid sections with the 
same tangents. Therefore the ellipsoid sections cannot be geodesic surfaces (in each 
point they have common tangents with different geodesic surfaces) . 
In order to have the generic geodesic motions 

d^a, , daj dak 

we need the Christoffel symbols r*j, of the space of orbits. They are given as 

P * P,9 

5^kaJ ^ ^liJpAipJ {Epj - Ejp) + {j ^ k) \ . (72) 
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In the above formula and the forthcoming ones all summations are indicated explicitly. 
The coefficients Ay^ and have the expressions; 

Ajjfe = al - 2af^ ruj + 2a]mk (73) 

r 

Y,al-2a^y^m, (74) 

r 

We note here that in the {ai} coordinates r*^, diverges on the boundaries. The geodesic 
equations (|7l|), supplemented by the boundary conditions (3) and (4) represent the 
equations of motion for the three-body problem with vanishing constants of motion, 
constant potential case, in terms of relative separations. 

We seek for solutions of the geodesic equation in the following form f 



X^t, a,{t) = y/A,t^ + B,t + a ■ (75) 

These characterize the uniform velocity, straight line motions of the particles in the 
physical space. For such motions the vanishing of the energy is assured by a proper 
choice of the constant V. The coefficients Ai, Bi and Ci are functions of the relative 
positions Yjk = Vj — and relative velocities Yjk of the three particles at some initial 
time: 



— 2 '^^l^ijk^'jk I — l^ijkl^jk '^jk , ~ 2 '^^l^^i'^^^ " ^"^^^ 

j.k j.k j,k 



Therefore Ai — and Ci — both imply Bi = too. 

It is straightforward to check that particular cases of the motions ( |75| ) fulfill 
the geodesic equation (|7l|). In the case Ai — the geodesic equations are trivially 
satisfied, this choice of the constants corresponding to no motion at all in the space 
of orbits. When d = 0, then ai{t) = \J~Ait are geodesies passing through the 
origin, with affine parameter i. These are the motions with tangents K"^ , the affine 
parameter t being linearly related to the radial coordinate r. The same type of motions 
ai(0 = -sTAit + \fCi emerge also for the choice Bi = 2\J AiCi. In this latter case Eqs. 
( [76| ) imply Ai = tpCi therefore Yjk = <prjfe- The special case of a pure expansion 
motion (a similarity transformation of the triangle of Fig. 2) is included here as the 
particular case with = ipVi. 

Those motions (^5|) which do not violate the constraints P = and L = are 
expected to be geodesies with affine parameter t. The fact that not all motions ( [75| ) 
are geodesies is readily seen for the case Bi ~ with all other constants nonvanishing. 
In this case the change of the relative velocities tjk is orthogonal to the relative 
positions Vjk ■ When all Ai are equal, this situation corresponds to an overall rotation 
of the system. 

The integral form of the generic geodesies is found as follows. Considering a motion 
of the type (|75| ) through a point ai{0) with direction ai{0) at t = 0, we find that the 
constants Bi and C'i are determined directly by the initial data 

B, = 2a,(0)d,(0) , a = a, (0)2 . (77) 



By enforcing the motions ( |75| ) to be geodesic through Eq. (|71|), three polynomial 
identities in t are found, each of degree 6. As the geodesic condition holds at all 

t For generic potentials V we do not dispose of such a solving Ansatz, as is well known from the 
three-body problem |7j. 
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instants, we can impose that all coefficients vanish. It turns out that the coefficients 
of t'^ contain the constants Ai linearly. Therefore, (or equivalently from the geodesic 
condition at t = 0) the constants Ai can be expressed in terms of Bi and Ci. It is 
then verified than with these Ai all other coefficients are also zero. The result of the 
above computation is 



Ai^l -mfCi (y^EijkAjjkBk) + 2mfBi (^^ei-jkAijkBk^ (^Cmi^n/ 

V j^k j^k n,l 



I j.k 



/4 (7(1)) ' {Cl + Cl + Cl - 2C1C2 - 2C2C3 - 2C3C1) , (78) 
where h.ijk and Sj^ are given by Eqs. (73), (74) and 7^") by 

7(") = dm^mj, + C^m-imi + Clmim^ . (79) 

The quantity 7^^^ is just 2MW introduced before. By taking the common denominator 
of Eq. (78), we find sums of expressions of the type AiCj and B^Cf on the left and 
right hand sides, respectively. It is easy to check, that all previously discussed special 
cases apply. 

The expressions (78) are indeterminate on the boundaries. Indeed, it can be shown 
that there the numerators are proportional to a^. The presence of the vanishing 
factor in the denominator is related to the singular behavior of r* j. in the coUinear 
configurations. 

Nevertheless, by starting the motion from any noncoUinear configuration, the 
coefficients Ai are well deffiied, and Eq. (75) together with the coefficients ( |77| ) and 



(78) represent the general solution of the geodesic equation in the constant potential 
case. With this, the Cauchy problem is also solved: we have found the motions 
pertinent to arbitrary noncollinear initial data (0^,0^). 



7. Concluding Remarks 

We have presented a detailed analysis of the geometry of the space of orbits of the 
Leibniz group for three particles. We have shown that in the line configurations 
(corresponding to the boundaries of the reduced space) no curvature singularity occurs. 
The boundary conditions (3) and (4) imply that the geodesies reach the boundaries 
tangentially after which they must return to the inner region. 

Our geometrical approach provides an alternative to the classical treatments of the 
three-body problem 0. As has been shown first by Lagrange, a generic reduction of 
the 18*'* order system of differential equations characterizing the three-body problem 
to a 6*'' order system is guaranteed by the existence of the ten integrals of motion. For 
generic situations no other reduction can be made, as the theorem of Bruns forbids 
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the existence of any further (algebraically independent) integral of motion. The 6*'' 
order system characterizing the particular motions with vanishing first integrals is the 
set of geodesic equations ([7l|). Supplemented by the boundary conditions (3) and 
(4), they describe free motions. We have investigated the constant potential case, 
however generalization is straightforward by a conformal rescaling of the metric, with 
the conformal factor —2V. We hope that the geometric methods developed in this 
paper will be applied in the study of nontrivial dynamical problems either, pertinent 
to specific V ^ const cases. 

For completeness we compute the Coriolis tensor for three particles ^ which was 
introduced initially in the context of molecular dynamics. In a body frame differing 
from the principal axis frame by a rotation about the third axis and choosing the basis 
dai in the shape space it has the components: 

Here the index v refers to the body frame and indices a, b and c to the shape space 
(our space of distances) . We would like to stress that the Coriolis tensor ( ^ ) is related 
in a simple way to the conformal Killing vector K"^, the existence of which we have 
exploited in many ways in this paper. 
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